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DEFINITION OF SYMBOLS 

Symbol 

l 

S/2 

Yield point in Shear 

o'r- o" 0 > cr z 

Radial, circumferential, axial stresses 

e r > e 0 > e z 

Radial, circumferential, axial strains 

y 

A function used in Hencky's theory 

e 

Y/H 

3H 

Bulk modulus of elasticity (corresponding to 3 a 
in Ref. 4 

E 

Young* s modulus of elasticity 

G 

Modulus of rigidity 

V 

Poisson's ratio 

m 

A numerical constant denoting the exponent of a 
simple power function 

r 

(variable) radius at a point of the cylinder 

r o 

Inner radius of the cylinder 

kr o 

Outer radius of the cylinder 

nr o 

Radius of elastic -plastic interface 

w 

r/r c 

1 f 

A function used in the modified criterion of yielding 

A, B, C 

Constants of integration 

Other symbols will be defined as they appear in the text. 
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STRESS AND STRAIN DISTRIBUTIONS IN A 
THICK-WALLED CYLINDER OF STRAIN- HARDENING MATERIAL 
ELASTIC -PLASTICALLY STRAINED BY INTERNAL PRESSURE 

By C. K. Liu* 

George C. Marshall Space Flight Center 


SUMMARY 


311*4 


An attempt is made to extend Allen and Sopwith's solution for an 
elastic -plastic thick-walled cylinder of ideally plastic material to 
strain-hardening material. Tresca's criterion of yielding is modified 
so that the effective stress (usually considered a function of the 
effective strain) is proportional to a simple power function of the radius 
in the plastic zone. An analysis of Allen and Sopwith's solution for 
ideally plastic material seems to support and justify the proposed 
modification of the yielding criterion. The exponent of the power 
function is taken as -3/2, -5/4, -1, and -1/2, in each case yielding 
an analytic solution corresponding to a particular strain-hardening 
material. Examples worked are for the case of plane strain and for 
a ratio of outer to inner radii of 2.0. When the exponent is taken to 
be zero, the solution becomes identical to that of Allen and Sopwith. 

For exponents other than the five mentioned above, numerical or ^ 
graphical methods of integration may be used. 



INTRODUCTION 

During the past two decades many solutions have been presented 
for the problem of plastically strained thick-walled cylinders. Because 
of different philosophies and purposes of the originators, some solutions 
emphasize the practical aspects of the problem (REF. 1, 2, 3); others 
are elegant in their mathematical operations (REF. 4, 5, 6); some are 
analytical (REF. 5, 6), and others numerical (REF. 7, 8, 9). 

F or engineering applications and for mathematical expediency, 
the authors impose various end conditions on the cylinder, or they 
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hypothesize that the material is incompressible. In some cases, where 
plastic yielding is predominant, the material is considered as being 
plastic-rigid. Moreover, the solutions also differ as the material is 
assumed to strain-harden or to be ideally plastic. 

Steele (REF. 6) presented a solution to this problem in closed and 
simple form. It is based on Hencky's total plastic strain theory, 
Tresca's yielding criterion, and the assumption that the material is 
incompressible in the elastic state. In assuming this, Steele postulates 
that across the entire thickness of the cylinder (in both the elastic and 
plastic zones) the strain decreases in inverse proportion to a function 
of the square of the radius at which the strain is sought. 

Earlier, Allen and Sopwith (REF. 4) effected an analytical solution 
to the problem. Since they include the compressibility of the material 
in the elastic zone, their solution is more rigorous mathematically, 
physically sound, and reasonably convenient to apply. Yet, it does 
not have provisions to take into account the strain-hardening property 
of the material. The solution in this report extends Allen and Sopwith 1 s 
solution to strain-hardening material by using a modified Tresca 
criterion of yielding. 

As shown by Hill, Lee and Tupper (REF. 7), the more exact theory 
governing the stress-strain relationship in the plastic region is the 
so-called "incremental theory" of which Hencky's total strain theory 
is a special case. The latter however, is the form used in the vast 
majority of work on plasticity (REF. 4). Furthermore, the solution 
obtained by applying Hencky's theory is a close approximation to that 
obtained by the more exact'theory (REF. 9). It is probable that these 
authors (REF. 2, 3, 5, 9), as well as Allen and Sopwith, have used 
Hencky's theory for no other reason than mathematical simplicity. 


METHOD OF SOLUTION 


The compressibility of the material can be expressed as follows: 


£ r + £ 0 + £ z, 3 (°r + 0-0 + °- z ) 


(1) 


where 


H = 


1 - 2 v 
3 E 


2 


By Hencky's total strain theory, the following relationship between 
the plastic strain and stress is postulated: 


6 9 ill. - € r - £ z _ e z - _ 3 y 

°e -<r r °-r - <r z “ <tq 


( 2 ) 


Tresca's yielding criterion is modified here as follows: 


°e - °r = S f 


( 3 ) 


where f is a certain function dependent on the strain-hardening 
property of the material. This function (f) must assume a value of 
unity when the strain-hardening effect of the material is absent, or 
when plastic yielding is impending. 


Continuity at elastic -plastic interface, in addition to the given 
conditions at the inner and outer cylinder surfaces, provides the 
necessary equations from which the stresses may be determined. 
These conditions are: 


cr r | elastic 


- cr r J plastic 


at 

w 

— 

n 

o-e ^r = 

S 

at 

w 


n 

<r T = 

O 

at 

w 


k 

<r r = 

-P 

at 

w 
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( 4 ) 


The stress distribution in the elastic zone is given by the Lame 
solution, that is, 


<r 


r 


°e 



A 



(5) 


W 


( 6 ) 


In the plastic zone ( 1 ^ w ^ n) } the stresses are found from the 
equation of equilibrium 


— CTq — CTj. — S f 


°r = S f — 
J W 


dw 4 - C 


= s 


[ f + fi dw ] 


+ c 


Applying boundary conditions (4) the constants A, B, and C are 
found to be 


A - — — 

2 k 2 


B = yn' 


C = S 


,2 _ 1,2 


) - 


Hence, the stresses across the thickness of the cylinder (FIG 1), are 
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FIGURE 1. AN INFINITE HOLLOW CYLINDER ELASTIC -PLASTICALLY 
STRAINED BY INTERNAL PRESSURE 



tr r n 2 1 1 “ I 
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elastic zone (n g w g k) (13) 


cr z E € z n- 

+ v r 

S S k 


/ w f 

— dw + n 2 - k 2 

w -g? — 


o-fi T f n 

-H- = f + / — + ~ 

S w 


2 - k 2 
2k 2 


-^2- = to be determined later 
S 


plastic zone (1 g w g n) (14) 


The internal pressure (p) can be found from condition (4) as 


r l {_ 
X w 


k 2 - n 2 


From the condition of compressibility (1) and Hencky's total strain 
theory (2), the strain in the plastic zone can be expressed in terms of 
the function Y and the stresses given by equations (13) and (14). 


(H + 2 V) tr z 

•f 

(H - Y) (o- r 

+ CTq) 

(16) 

(H + 2 Y)<r r 

+ 

(H - Y) (< t z 

+ <r 0 ) 

(17) 

(H + 2 Y )o* e 


(H - Y) (<r r 

+ ^z) 

(18) 
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Rewriting, equation (16) becomes 


0-7. = 


H + 2Y 


H - V 

H + 2Y • + °e) 


( 19 ) 


Substitution of equation (14) into (19) gives 


= H+ Z 2Y • jfrw [ f + z j[ £- dw + ] 

Consequently, the three components of strain are 

| ~3Y(2H + Y) f 
[ H + 2 Y 

1HY /r W f , n 2 - k 2 M 

f 3Y(H -Y) , 

[ H + 2Y 

9hy ( r w x , , n 2 - k 2 Vt 
+ 2Y \J n w dw + 2k 3 /J 


2 - k 2 

5— | S (20) 


€ e = h + 2Y € * + s 

_9HY 
+ H 

H - Y 

€r = H + 2Y €z + S 


e z = constant, to be determined from the end condition of the 
cylinder 


( 21 ) 


( 22 ) 


(23) 


To further investigate the relationship among the strains and the 
function f, it is necessary to apply the equation of compatibility 


w 


d € fl 
d w 


- € 


e 


(24) 


The following equation is obtained after equations (21) and (22) are 
substituted into equation (24) 


/ dw 

df \ 

dY 

2 V 

+ t) 

= Y 


_f 


H € r 


- 3H 


CT- 

Wn w 


dw + 


n 2 - k 


2k 2 




2(H + HY + Y ) (25) 
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Now, if f = 1 in equation (25) and we write 


Y 



-j- = 2 ln(w/n) 


n“ 

F 


2 

3 SH 


( 26 ) 


it is found that equation (25) reduces to 

dy 3y _ _ 3x 

dx + x (x - 1) (x + 3) (x - 1) (x + 3) 


(27) 


Equation (27) is identical to that of Allen and Sopwith (REF. 4, p. 70) 
for an ideally plastic material. Its solution is found as 


ln(w/ n) = 


3 HS 


m 

2P 


1 + 3u 4 
iu 3 


“ /3n 2 

“ ( "F " 


2 

3 




For convenience in later computations, it is assumed that e z = 0 
(that is, the cylinder has clamped ends) and that v = 1/3. Thus, 

equation (28) is simplified greatly: 
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n x 


ln(w/n) = - + 


1 - 3u 4 


- tanh”’(2/ 3) 


lu 5 

l/4l 



7/4 


t-l 


tan -1 (2/3) 


1/4 


+ (tan -1 u - tanh‘ 


" u, l 


( 29 ) 


Equation (29) correlates the quantities w/n and u (or Y , in turn) 
in a transcendental manner. The term n/k appears to be a parameter 
that may or may not affect the (w/n) versus Y relationship appreciably. 

To determine this, w/n is plotted against Y/H (= £) in FIG 2 for a 
series of values of n/k: 1.0, 0.875, 0.750, 0.625, 0.500. In the 
actual cylinder, w lies between 1 and n in the plastic zone. The 
corresponding range of w/n is 1/n to 1. Hence the plot of Y versus 
w/n for a particular value of n/k has physical significance only when 
w/n is within the range [l / k(n/k)] to 1. This is the way to determine 
the range of interest of w/n for any choice of k. 

In this report the case of k = 2 is considered. It is found that 
while n/k varies from 1 to 0. 5, the maximum range of variation in 
w/n never exceeds 0. 37 percent. With the scale of plotting used in 
FlG 2, it is difficult to distinguish these four curves, one overlapping 
whose relationship with Y is not appreciably affected by a change in the 
value of n for a given value of k. 

Considering w/n as an independent variable, it may be advantageous 
to assume that the function f is a function of w/n only. A possible 
form of f may be written as 

f = C,(w /n) m + C 2 (w/n) a + C 3 (w/n)^ 

+ ( 1 - w ^ n) (30) 

where C } , C 2 , C 3 , . . . . ; m, a, p, . . . . are constants. The purpose 

of shaping f in this manner is to simulate the unknown relationship 
between w/n and the respective plastic strain. Here, according to 
the characteristics of f, it will be necessary that when w/n = 1 

C 1 + C 2 + c 3 + =1 (31) 


Q 
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FIGURE 2. A DIMENSIONLESS PLOT SHOWING RELATIONSHIP 
BETWEEN £ AND — FOR VARIOUS VALUES OF 

VI Lr 




Let it be further assumed that the constant C* is predominantly 
larger than the sum of the rest of the terms, then Cj = 1, 

C 2 = C j — # • • • — 0 and 


f a (w/n) m 


(32) 


With the function f expressed by equation (32), equation (25) 
becomes 


df (3 + 2m) H 2 + 2mHY 4- 2m V 2 
dY + (2 + m) Y (H + 2Y) (2H + 2Y) 


m H j 6z 

m + 2 Y(2H + Y) (H + 2Y) |_ S 




(33) 


The problem is thus reduced to solving the first order, ordinary 
differential equation (33). The solution of equation (33) is 


/w T 1 1 + 2§ | 

m j£ z 

\n) * 2m + 3 


2m + 1 

m + 2 |_SH 

^2(m + 2) 


2(m + 2) 


(2 + i) 



3 /m + 2 

n 2\l 



+ 2 \ m 

f 

■ k^/J 

dt 



1 


+ K' 


2(m + 2) 


e 


( 1+2 £ ) 2 (2 + £) 


2(m + 2) 


( 34 ) 


in which £ = Y/H, and K', the constant of integration, is determined 

by the condition: w = n, or f = 1, £ = l/(6GH). Thus: 

2m + 3 Zm + 1 


/ i \ 2(M + 2) / j \ 

K , „ Ugh) V^gh/- 


2(m + 2) 


m 


m + 2 


1 + 


3GH 


Ja. + 1 (™ + 2 n 2 \ 

SH 2 \ m ky_ 



de 


i 


2(m + 2) 


2(m + 2) 


(1 + 2 £) 2 (2 +|) 


(35) 


1 

£ = 6GH 


Consequently, equation (34) becomes 


1 + 24 
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2m + 1 




2(m + 2) 

m _ 5 (2 + 4 ) 

2(m + 2) 

m 


,3 /m + 2 n 2 

W 

i i A 


m + 2 

. SH 

+ 2 \ m k 7 

3GH 


§ 


1+24 

2m + 3 

2m + 1 

2m 

+ 3 

2m + 1 

2(m + 2) 

/ 1 1 L . 

2(m + 2) 
1 \ 

e 2<m 

+ 2) 

2(m + 2) 

(2 +4) 

\bGH/ 1 6GH J 



r 5 

d 4’ 





3 


+ 2) 2(m + 2) 

i * (i + 2 e 1 ) 2 (2 + e') 

6GH 



(36) 
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The evaluation of the integral in equation (36) for the general case 
is quite cumbersome, but the integration can be carried out easily if 
the exponent m is given particular values, such as -3/2, -5/4, -1, 
-1/2, etc. Four solutions corresponding to these four exponents are 
given below, 
m = -3/2, 


U + 2j\ (2 +S) 2 

(* + 3Gh) ( 2 + 6GH/ 

j' + 22 ., 1 _ + 

g' + 2 27 (2 + £') 




2 1_ 

9 (2 + fj 7 





(1 + 2 1 ) (2 +£) 2 . 


1 

(1 + 2 S') 



6GH 


(37) 


m = -5/4, 



(1 + 2g) (2 + £) 

,1/3 



1 i£_ 

3 [SH 


t nl/3 1/3 

4 e 1 §' 

9 1 + 2 £' + 18 2+1' 




+ 


+ 


( 2) 1/3 + €' 1 /3 

7T 


27(2) 2/3 ln C.2/3. 1/3 . 1/3 2/3 11/2 

' ' t* -(2) $ +(2) ; 


5(3) 


1/2 


27 (2) 


2/3 


1/3 

tan"’ 2JLT_M2) 


+( 2 ) 

1/3" 


(3) 1/3 (2) 1/3 


(38) 


6GH 



m = - 1, 





1 / e 1 - 2 

9 v! 73 - F + 1 



1 + S +[£ (g + 2)] 1/2 

1 + 6GH + (1 + 12GH) 1/2 (39) 

~ ~ 6GH 



l_3 1 + 2 g' 9(2) 


1/3 


In 


le 2/3 1/2 

( 2) 1 /3 + £• 1/3 


(3) 

9(2) 


1/2 

ITT 


tan 


-i 2 g' 


1/3 


“ ( 2 ) 


1/3 




(2) 1/3 (3) 1/2 


6GH 


(40) 
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Again, assuming « z = 0 for a cylinder with clamped ends, and 
v = 1/3 (6GH = 1/4) it is possible to examine the relationship between 
w/n and £ for various values of m . In FIG 3, four curves are 
plotted so that the effect of m on £ can be brought out. It is seen 
that the material with m = -3/2 at a fixed value of w/n yields the 
smallest value of £ among the four materials. Note that fixed value 
of n corresponds to a fixed magnitude of the internal pressure (p) 
through equation (15). In FIG 2, a dimensionless plot is shown for the 
relationship between £ and w/n for various values of n/k. The 
single curve shown is actually four curves of different lengths over- 
lapping one another. 

With the interdependence between (w/n) m and £ determined, it 
is possible to check back to see what degree of strain-hardening each 
value of m represents. From the yielding criterion (3) it is assumed 

Ffl Z = f (3) 

S 


and from Hencky’s total strain theory (2), 


6 e - 6 r 
- £ r r 


3Y 


( 2 ) 


Combining equations (3) and (2), 


Ifi L_ 

S H 


€ r 



(41) 


For certain values of m, £ and f are related to each other 
according to the curves in FIG 3. Hence, by plotting f against 3£f, 
a dimensionless stress-strain diagram of this particular material can 
be obtained. Four curves of this kind, including that for an ideally 
plastic material (m = 0, are shown in FIG 4. 
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* Once the exponent of strain-hardening of a cylinder material is 

selected by comparing its dimensionless stress - strain diagram with 
I the curves in FIG 4, the stress and strain distributions across the 

cylinder radius can be readily determined from equations (13), (14), 
(19), (21) and (22). The stress distributions in cylinders made of four 
different materials (m = 0, -1/2, -1, -3/2) for a degree of yielding 
of n = 2, are plotted in FIG 5. The cylinder used in the illustrative 
example has an outside radius twice as large as its inner radius, in 
other words, k = 2. 

A comparison between equation (19) and the third of equations (13) 
shows that the factor 

V - H £ - 1 

H + 2Y ° r 1 + 2 £ 


has characteristics similar to that of the Poisson's ratio, v. This 
factor assumes a value of 1/3 when £ = 4, and approaches the value 
of 1/2 when £, — ► go. 

CO NC LUSIONS 

Allen and Sopwith's solution for the thick-walled cylinder problem 
for an ideally plastic material was extended to strain-hardening 
material by introducing a strain-hardening function to the Tresca's 
criterion of plastic yielding. This strain-hardening function is assumed 
to have the form of a simple power function of the dimensionless radius 
in the plastic zone. The exponent of this simple power function may 
have different constant values, fractional or integral, depending on the 
strain-hardening property of the material concerned. As demonstrated 
in the foregoing analysis, the solution to this problem is in analytical 
and closed form if this exponential constant is equal to -3/2, -5/4, -1, 
-1/2 and 0. These exponents correspond to five materials with different 
strain-hardening characteristics similar to those observed in metals. 
Since the case of zero exponent coincides with that of Allen and Sopwith's 
ideally plastic material, their relative degree of strain-hardening can be 
compared by examining the curves in FIG 4, where the dimensionless 
quantities of maximum shearing stress and strain are plotted as 
ordinate and abscissa. 
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Stress/ S 


n = 2 



FIGURE 5. FULLY PLASTIC STRESS DISTRIBUTION IN AN INFINITE 
HOLLOW CYLINDER STRAINED BY INTERNAL PRESSURE 
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The assumption that the maximum shearing stress (effective stress) 
is proportional to a simple power function of the radius has been deduced 
from the solution of Allen and Sopwith as reasonable for engineering 
applications. Examination reveals that the value of this simple power 
function at any radius in the plastic zone is not appreciably affected by 
the amount of yielding that has penetrated beyond that redius. Within 
the realm of engineering significance and application, the example 
shown in this report (k = 2) seems to strengthen this assumption. 

For exponents other than the five mentioned above, the integration 
involved may become very tedious for analytical solutions. In those 
cases, graphical or numerical processes may be used. 

Although the cylinder used in the illustrative example was assumed 
to have clamped ends, this was mainly for the convenience of computation 
and is not a limitation of the theory. 
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